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EFFECT OF STATIC STRESS AND EDGE RESTRAINT ON THE VIBRATION 
OF NEARLY CYLINDRICAL SHELLS WITH CIRCULAR CROSS SECTION 

By Paul A. Cooper and Martha P. Robinson 
Langley Research Center 

SUMMARY 

The effect of stress due to statically applied axial compressive and tensile loads 
and external or internal pressure on the minimum natural frequencies of freely supported 
nearly cylindrical shells with circular cross section and slight meridional curvature 
is investigated. Also the effect of various edge restraints on the natural vibration of 
unstressed nearly cylindrical shells is investigated. Results show that axial tension loads 
and external hydrostatic pressure can cause sizable reductions in minimum frequencies 
for positive Gaussian curvature shells. While axial compression loads and internal pres- 
sure increase the frequencies for these shells, the compression loading increases fre- 
quencies only slightly before onset of static instability. In contrast, negative -curvature 
shells are highly sensitive to axial compression and external hydrostatic pressure and 
evidence large reductions in frequencies at low load levels. Both axial tension and inter- 
nal hydrostatic pressure cause large increases in frequencies in negative -curvature 
shells, but generally not higher than similarly loaded cylinders. 

The minimum frequencies of negative -curvature shells are sensitive to the degree 
of meridional in-plane restraint at the shell edges, with large increases occurring when 
the edges are completely restrained, whereas they are relatively insensitive to slope 
restraint at the edges. The minimum frequencies of positive -curvature shells, on the 
other hand, are relatively insensitive to both in -plane and slope restraints. 

INTRODUCTION 

Doubly curved shell elements are commonly used in aerospace vehicle structures. 
With the exception of the special case of spherical shells, very little general information 
on the dynamic behavior of doubly curved shells exists in the literature. One class of 
shells which is useful for characterizing the behavior of doubly curved shells is the nearly 
cylindrical shell with circular cross section which contains slight meridional curvature. 
The vibration behavior of such shells has been studied analytically in reference 1, for 
example, in which unstressed shells with "freely supported" edges are treated. The term 
freely supported denotes edges which have zero displacement in the direction normal to 
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the shell middle surface but which are free to displace in the direction along the middle 
surface meridian and free to rotate. It was found in reference 1 that for nearly cylin- 
drical freely supported shells with negative Gaussian curvature, the minimum natural 
frequencies are generally below those of the corresponding cylinder and evidence wide 
variations in value as the meridional curvature is varied. At certain critical values of 
negative curvature, exceptionally large reductions in frequency occur. These large 
reductions result from a condition of essentially inextensional deformation during the 
vibration. Positive curvature shells, on the other hand, exhibit higher minimum frequen- 
cies than those of the corresponding cylinder. 

The purpose of the present paper is to extend the investigation of reference 1 in 
two ways. First, the effects of several constant biaxial membrane stress fields on the 
natural vibration of freely supported shells are presented. Second, the effects of various 
types of edge support as means of preventing inextensional vibration modes are investi- 
gated for unstressed shells. As in reference 1, the minimum natural frequency is used 
as an indicator of the vibration behavior of the doubly curved shells. 

The governing equations used for this analysis were developed in reference 1 and 
are restricted to nearly cylindrical thin shells having shallow meridional curvature. The 
equations are solved for a general linear elastic edge restraint by methods similar to 
those reported in references 2 and 3, where trial values of frequency are introduced to 
force the vanishing of a characteristic determinant. The method of solution is "exact" in 
the sense that the natural frequencies can be found to any degree of accuracy desired. 

SYMBOLS 

Aj displacement coefficients defined by equations (5) 

B extensional stiffness 

c central rise of the shell meridian 

D bending stiffness 

E Young’s modulus of elasticity 

hj,gj displacement coefficients (eqs. (7)) 

h shell thickness 

m number of axial half-waves 
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M|,M 0 ,M £0 moment resultants associated with vibration state 

n number of circumferential waves 

N£,N 0 ,N £0 stress resultants associated with vibration state 

N£,Nfl static stress resultants 

p,N surface loading, edge loading 

r circumferential radius of doubly curved shell 

R radius of cylinder; circumferential radius at midlength of doubly curved shell 

meridional radius of curvature 
s total meridional arc length 

t time 

u,v,w displacement variables in meridional (£), circumferential (6), and normal 

directions, respectively, defining the vibration state 

U,V,W coefficients in equations ( 8 ) 

a,fi,y,5 coefficients of characteristic roots (eqs. ( 10 )) 

8 circumferential coordinate 

Xj characteristic roots 

X thickness parameter, h/R 

Aj coefficients of equation (20) 

p Poisson’s ratio (taken as p = 0.3 for all calculations) 

v mass density 
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Nfc . 

nondimensional meridional stress, x 10’ 


meridional coordinate 


percent ratio of meridional rise to length, §(100) * ^—(100) 

& oKfc 


0) 

0 

V4 

Matrices: 


natural circular frequency 


frequency parameter, Ru>| 


Ml - M 2 ) 

E 


A q4 n q 4 g4 

biharmonic differential operator R* — -r + 2R ,i — £ — « + — t 


n- 


9 £ 2 90 2 9 #4 


{b) 8x 1 colvunn matrix 

&],[n],[k],[y] 8x8 square matrices 

Subscripts following commas represent derivatives with respect to variables indi- 
cated by the subscript. 


DESCRIPTION OF ANALYSIS 


Governing Equations 

The equations governing the motion of symmetrically stressed nearly cylindrical 
thin elastic shells of revolution with shallow meridional constant curvature of the type 
shown in figure 1 are given in terms of displacements in reference 1 (with origin of | 
taken at midlength) as 
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Negative Gaussian curvature shell 


Figure 1.- Geometry of class of shell investigated. 
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where 


X = h 


R 


and the circumferential radius is defined by r = R - 


2R 


i 


The quantities u, v, and w 


are the displacements in the meridional, circumferential, and normal directions of the 
shell surface, respectively; R is the circumferential radius of the shell at midlength; 
R^ is the constant meridional radius of curvature; h is the shell thickness; and \i is 
Poisson's ratio. The variable £ is the measure of length along the meridian with its 
origin taken at midlength (fig. 1), 6 is the angular circumferential measure, and t is 

time. The quantity v is the mass density of the material and B is the extensional 
stiffness of the shell surface given by 


R — Eh 


( 2 ) 


where E is Young's modulus of elasticity. Equations (1) have been derived on the 
assumption that the shell material is homogeneous, isotropic, and linearly elastic. 
Transverse shear and rotary inertia terms have been neglected. 


Deformations due to initial static loading have been assumed to be uniform, and 
variation in stress along the meridian due to double curvature is zero in the approximate 
theory, so that the static stress quantities and Ng are constant membrane stress 
resultants given by 




= N 



(3) 


where N is an applied tensile edge meridional force resultant and p is an applied 
uniform internal constant directional pressure. 


The corresponding boundary conditions at £ = are given by 




(i £♦§)-» 
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■ W = 0 


3 3 , 91 ; 

s [-< 2 - “> v ,{e - r2w ,|{{ + < 2 - "> w ,eeJ - 


N & 

i^.r 0 




or u = 0 


or v = 0 


or w = 0 


or w t = 0 1 
> s 


( 4 ) 
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Stressed Shells With General Boundary Conditions 

Equations (1) are constant-coefficient partial -differential equations which are sat- 
isfied by a solution of the form 


D 

u = Ajhje sin n0 cos art 

D 

v = Ajgje cos n 6 cos art 
w = Aje R sin n0 cos art 




(5) 


where n is the number of circumferential waves, a > is the natural frequency, and the 
Xj's are the characteristic roots of equations (1) and must satisfy the biquartic charac- 
teristic equation 


a 6 Xj 6 + a 4 Xj 4 + a 2 Xj 2 + a 0 = 0 


( 6 ) 


The coefficients ag, a 2 , a 4 , and ag are given in appendix A and are functions of p, 
n, R, Ft£, and £2. 

The amplitude coefficients hj and gj of equations (5) may be found from any two 
of equations (1). From equations (lb) and (lc), they become 




■( f 6 x j 2 + *7) + ( f 8 X j 4 + f 9 X j 2 + f lo) ( f 4 x j 2 + *5) 

4 



•f2(f8 x j 4 

+ fgXj 2 + flo) + f 3 6 X j 2 + f 7y 



f 2( f 6 x j 2 +f 7) “ f 3( f 4 x j 2 +f 5) 






J 


( 7 ) 


where the functions fj, . . ., f^g are S iven in appendix A. 

For the special case of axisymmetric vibrations (n = 0), the terms gj vanish and 
equation (lb), the circumferential equilibrium equation, uncouples from the remaining 
equations (1) so that the procedure must be modified. The axisymmetric case is handled 
separately in appendix B. 
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The complete solution for u, v, and w is found by summing the contributions 
from each of the eight characteristic roots of equation (6) so that 


u = U(£)sin nd cos <ut~' 
v = V(£)cos n 6 cos a>t / 
w = W(£)sin n0 cos cot, 


( 8 ) 


where 


8 

u(i) = 1 Ajhje X i 5 
i=i 

8 

V(£) = ^ A jgje X ^ / (9) 

j=l 

8 

W(£) = £ Aje^ 1 
3=1 


For the problems treated in this analysis, the natural frequencies of interest are 
always associated with characteristic roots of equation (6) of the following form: 

Case (1): 


Aj = ±(a ± i0), ±y, ±i6 

(10a) 

Aj = ±(at ± i/3), ±iy, ±i6 

(10b) 


where a, y, and 6 are real and in general distinct. The natural frequencies of a 
cylinder were generally found to occur for characteristic roots of the form of case (1) in 
reference 2. In the present analysis, case (2) can also occur. 

The displacement coefficients given in equations (9) can be written after some 
manipulation in terms of real functions as 
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r — \ 

U(£) = 

Z B jUj(l) 


j=l 


8 
1 \ 

V(£) = 

X 


3=1 

8 

W (|) = £ BjWjfe) 


(ID 


where Uj, Vj, and Wj are real functions of £ given in appendix A. The amplitude 
coefficients Bj are undetermined real constants which take the role of the eigenvector 
elements in the vibration solution. 

With the origin of £ taken at the midpoint of the shell, the boundary conditions of 
equations (4) may be written in general at two edges £ = ^ and £ = --| as 



where {b} is an eight -element column matrix for which the elements are Bj and 
where |n] and |V] are 8X8 square matrices for which the elements are given in 
appendix A and where [$] and [k] are 8x8 square boundary selection matrices for 
which the elements are defined as 


*jk = K jk = 0 


(for j * k) (13a) 


and 


K jk = 1 - % (for j = k) (13b) 

The elements 4>jj take on the value one or zero according to the prescribed edge condi- 
tions. If a stress -free condition is prescribed, the corresponding €>jj term is set equal 
to one. If a displacement variable or meridional slope is fully constrained, the corre- 
sponding <3>jj term is set equal to zero. The [¥] and [k] matrices may be general- 
ized to enforce coupled linear elastic or directional constraints. 
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For a given set of homogeneous boundary conditions and a particular circumferen- 
tial mode number n, the natural frequencies of the system are those which cause the 
determinant of the coefficient matrix in equation (12) to vanish; that is, 


HM + HH 


= 0 


(14) 


There are an infinite number of natural frequencies associated with each circumferential 
mode number, with each frequency in general associated with a unique natural mode. 
Because of the complicated transcendental character of the determinantal equation (14), 
a trial procedure is used to find the natural frequencies. The procedure is exact, how- 
ever, in the sense that the frequency parameter can be found to within any desired degree 
of accuracy. 

With a particular circumferential mode number n given and a trial value for the 
nondimensional frequency selected, the characteristic roots are determined from 
equation (6), hj and gj are found from equations (7) and Uj, Vj, and Wj are found 
from equations (A2) and (A3) in appendix A. With the proper case number determined, 
the elements of [n] and [y] needed to form the boundary conditions are calculated. 
The determinant of equation (14) is then evaluated and the result compared with zero. A 
new trial value is selected for Q, and the procedure is repeated until a change in sign 
of the determinant is observed. A Newton-Raphson convergence routine is then employed 
to select a new S2 and the same procedure is repeated until equation (14) is satisfied to 
within the desired accuracy. 

The success of this method as a rapid means of obtaining solutions depends on the 
selection of a good initial estimate of the frequency and a proper increment in frequency 
such that two roots are not traversed in one step. The closed-form solution for the nat- 
ural vibration of nearly cylindrical shells with freely supported edges (see next section) 
can be used to good advantage as an aid in the selection of initial trial frequencies for 
systems with other boundary conditions. 

When a natural frequency has been determined, its corresponding mode may be 
found from equations (11) and (12) . The 7x8 augmented matrix Cjj is formed by 
deleting one of the dependent equations from equation (12). With Aj defined to be the 
determinant formed from the square matrix found by deleting the jth column of Cjj, the 
modal amplitude coefficients Bj are proportional to (-l)j + *Aj and can thus be defined 
as 


B i = 


(-D 


j+1 


(15) 
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where e is some arbitrary scale factor. With Bj known, the meridional mode corre- 
sponding to S2 may be found from equations (11). Similarly, the vibratory stress and 
moment resultant distributions become, from the definitions of reference 1 and equa- 
tions (8) and (11): 


where 


D 

n £ = ft n £ sin ne cos 


•D 

Ng = g ng sin n0 cos cot 


LI 

N| q = — n^ g cos nd cos cot 


m^ sin n0 cos cot 
R“ 


Mg = - mg sin n 0 cos cot 
R z 

g = - m^ g cos nd cos cot 
R 


(16) 


and 


D = 


Eh 3 

12(l - ix 2 ) 


(17) 


n { “ 1 B i 
i=i 




8 


n e = 1 


J-i 

8 


M Uj' - nV. + (1 + n g-)W 


Re "3 


n ^= I B irin{ nls i + v i) 

3=1 

8 

mg = Y B jt nV j + w j" - Jun 2 Wj] 

3=1 ^ ' 


(18) 


Equations continued on next page 
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8 

£ Bj(-nVj + MWj" - » 2 Wj) 
i=l 
8 

m { „= £ Bj(l - WnWj' - Tj') 

3=1 J 

The primes on symbols in these equations indicate differentiation with respect to £. 


Stressed Shells With Freely Supported Edges 

A simplified procedure of solution is available for one set of boundary conditions. 
The exact solution of equations (1) for freely supported boundary conditions 
= v = w = = o) is given by 


u = U mn cos ^ ^ sin n 9 cos cot 


v = V mn sin 


'Ml., 

S 

hi) 


-\ 


cos n0 cos cot 


(19) 


’(« +!) . 

' ZL cir 


mn7 

w = W mn sin ^ — — sin n0 cos cot 

o 


where m is the number of meridional half waves, n is the number of circumferential 
waves, and co is the natural circular frequency. Substitution of equations (19) into equa- 
tions (1) leads to the following characteristic equation for natural frequencies: 


-SI 6 + A 2 S2 4 - A X S2 2 + A 0 = 0 


( 20 ) 


where S2 is a nondimensional frequency given by 


0 = 



( 21 ) 


The coefficients Aj depend upon the shell geometry, static stresses, and mode shape 
and are given in appendix A. 
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RESULTS AND DISCUSSION 

The analysis just described has been used to investigate effects of various stress 
states and edge restraints on the vibration behavior of doubly curved shells. Stress 
states treated include those caused by constant directional axial loads and internal and 
external hydrostatic pressure on vibrating freely supported shells. Edge restraints 
treated include zero in-plane displacement and zero slope on unstressed vibrating shells. 
All results presented are for nearly cylindrical shells having both positive and negative 
Gaussian curvature with a ratio of circumferential radius to thickness of 500, ratio of 
length to circumferential radius of 2, and Poisson’s ratio of 0.3. 

Results are presented in terms of an auxiliary parameter t, the ratio of the merid- 
ional rise to shell length expressed in percent: 

r = §(100) 

where c is the central rise of the shell meridian. (See fig. 1.) For small percentages 
of central rise with respect to length, r can be closely approximated by 

r = gf^(lOO) (22) 

It is shown in reference 1 that for unstressed shells the approximate theory from which 
the basic equations (1) are derived gives reasonable estimates of frequencies and predicts 
trends adequately for values of r between ±5 percent. 

Freely Supported Stressed Shells 

Axial tension . - Effects of axial tension on the vibration of freely supported shells 
of double curvature are shown in figures 2 and 3. In figure 2, the envelope of frequencies 
are plotted as a function of the ratio of meridional rise to shell length r, for various val- 
ues of an in-plane stress parameter a, which is essentially the meridional stress resul- 
tant normalized with respect to the in-plane axial stiffness of the shell wall (a = ^ x 10^. 

As the positive curvature increases from zero (cylinder) on the frequency envelope 
for cr = 0, the lowest natural frequency and its associated circumferential wave number 
increase. In reference 1 it was shown that the membrane frequency (where the shell is 
regarded as a membrane) can vanish for unstressed freely supported shells when 
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Minimum nondimensional frequency. 




Minimum none 



Figure 3.- Effect of axial tension on the minimum frequencies of several freely supported nearly cylindrical 
shells. The meridional wave number is 1; the circumferential wave number n is given on curves. 



Thus, unlike positive -curvature shells and finite-length cylinders, negative -curvature 
shells can develop an inextensional vibrational mode. As the negative curvature increases 
from zero, the minimum frequency decreases and occurs in a progressively lower cir- 
cumferential mode. For each circumferential mode, there is a distinct minimum in the 
variation of minimum frequency with curvature. These minimums occur at specific val- 
ues of curvature for which membrane theory predicts inextensional behavior in the asso- 
ciated mode. (See ref. 1.) The frequencies at these minimums are governed almost 
entirely by bending stiffness. 

Thus, for specific combinations of modal wavelength ratio and meridional curvature, 
negative -curvature shells vibrate in a nearly inextensional mode and have effective vibra- 
tion stiffnesses similar to those of simply supported plates vibrating with the same wave- 
length ratios. Reference 1 discusses in more detail the effect of curvature on the mini- 
mum natural frequencies for unstressed freely supported shells. The natural frequency 
envelope for a = 0 is repeated in all subsequent figures of this type for purposes of 
comparison. 

For an applied tensile stress resultant N at the shell edges, a constant biaxial 
membrane stress field is developed. From equations (3), the circumferential and merid- 
ional stress resultants are 


= N 




(24) 


The axial -tension load causes tensile stresses in both the meridional and circum- 
ferential directions in the negative -curvature shells which introduce a stiffening effect. 
As the loading increases, this biaxial tensile state eventually drives the minimum fre- 
quencies up above that of a similarly loaded cylinder (a > 11 in fig. 3). For the higher 
loadings, there is only a very slight variation in frequency with the degree of curvature, 
contrasted with the large variation that occurred with no load. 

For positive -curvature shells, on the other hand, the axial tensile edge load intro- 
duces compressive circumferential stresses which cause reductions in the minimum fre- 
quencies. Increased loading on these shells may lead to large enough circumferential 
compression to cause buckling. This condition is evidenced in figure 2 where the curves 
for a = 15 and 100 become almost vertical (i.e., S2 approaches zero) in the positive- 
curvature -shell regime. 

From the third of equations (1) and equations (24) it can be shown that the natural 
frequencies of a nearly cylindrical shell are independent of applied axial tensile or com- 
pressive edge loading when the associated mode has the relationship 
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(25) 




w ,00 

w M = mq 

if w = 0 at the edges. For a freely supported shell, this condition is satisfied when 
R _ 

R| \ ns / 

Every positive -curvature shell, therefore, has modes of vibration with a specific 
wavelength ratio for which the corresponding natural frequency is independent of the 
static stress. These modes of vibration are usually not associated with the minimum 
frequency of the shell. For the particular shell configuration treated in this report, 
however, this condition does occur on the minimum frequency envelope in figure 2 where 
the a - Constant curves intersect with m = 1 and n = 9 for r = 0.76 percent when 
ct is approximately equal to or less than 10. 

Figure 3 is a plot of the lowest natural frequency as a function of meridional non- 
dimensional stress for shells with various curvatures under axial tension. The minimum 
frequencies of positive curvature shells are very nearly insensitive to stress changes at 
low stress levels but become highly sensitive to changes in stress as the instability stress 
level is approached. 

The negative-curvature-shell frequencies are almost linear with stress on the log- 
log plot of figure 3 so that a simple power relationship (S2 cc a a ) would be a good indicator 
of the effect of axial tension. Figure 3 clearly displays curvature independence at a non- 
dimensional meridional stress level of approximately a = 11. At this stress level the 
natural frequency is nearly independent of the degree of curvature; hence, this level can 
be interpreted as the level where the stiffening effect of positive curvature is offset by 
the weakening effect of circumferential compressive stress. This effect can be observed 
in figure 2 where for a = 10 the minimum frequency envelope is almost horizontal for 
most of the range of curvature, the influence of curvature being felt only in the upper 
range of positive curvature. 

Axial compression. - Figure 4 is a set of envelopes of minimum natural frequencies 
plotted as a function of the central rise -length ratio for a range of nondimensional stati- 
cally applied edge compressive stress resultants a. The axial compressive load causes 
compressive meridional and circumferential stresses in the negative -curvature shells 
which severely weaken the shell and cause large reductions in frequency in local ranges 
of negative curvature and loss of stability at levels (ct = -0.5) an order of magnitude below 
the instability level for the corresponding cylinder (o' = -10.9). 
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In figure 5, the minimum natural frequency is plotted as a function of meridional 
compressive nondimensional stress for shells with various degrees of curvature. The 
axial compression introduces a static positive circumferential stress which tends to 
stiffen slightly the positive -curvature shells, but in general the effect is small. The 
positive -curvature shell becomes unstable at the same stress level as the cylinder 



Figure 5.- Effect of axial compression on the minimum frequencies of several freely sup- 
ported nearly cylindrical shells. The meridional wave number is 1} the circumferential 
wave number n is given on curves. 
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(ref. 4). The curves end abruptly at the buckling stress of the cylinder and do not exhibit 
a continuous decrease in frequency when the buckling stress is approached as does the 
cylinder. The minimum frequencies of the negative -curvature shells show more depen- 
dence on stress at low stress levels than do those of the positive -curvature shells and 
are highly dependent on stress only at stress levels close to instability. The particular 
curvatures at which membrane theory predicts inextensional behavior are shown as 
dashed lines and indicate the extreme sensitivity of shells to stress in this region of 
curvature. 


Internal hydrostati c pressure.- Figure 6 shows the effect of curvature on the mini- 
mum natural frequencies of internally pressurized shells. The internal pressure intro- 
duces a tensile meridional stress as well as circumferential stress in both positive- 
and negative -curvature shells and increases the stiffnesses of both shells above that of 


unstressed shells. For pressurized shells 


- 

a = 


B 


x 10*1 the envelope of natural fre- 


quencies in the negative-curvature range is consistently lower than that of the cylinder, 
although the difference in minimum frequency is no longer as large as for the unstressed 
shell. The positive -curvature range does not show large increases in stiffness. Thus 
the internal pressure acts as a stabilizing influence and generally moderates any large 
effects due to meridional curvature. 


External hydrostatic pressure . - Figure 7 shows the effect of curvature on the lowest 
natural frequency of externally pressurized nearly cylindrical shells. External pressure 
introduces compressive meridional and circumferential stresses. The effect of static 
external pressure on the lowest frequencies of negative -curvature shells is very similar 
to the effect of axial compression. Large reductions of frequency occur at curvatures 
corresponding to nearly inextensional vibration. 

Figure 8 shows the buckling pressure for negative -curvature shells to be well below 
that of the corresponding cylinder (on the order of 50 percent lower for r = -1). Unlike 
the results under axial compression, the positive -curvature shells are stable at pressures 
higher than the buckling pressure of the corresponding cylinder. The frequencies for both 
the positive- and negative -curvature shells are essentially insensitive to the pressure at 
low stress levels but are highly sensitive to stress as their respective buckling pressures 
are approached. 


Effects of Edge Restraint 

The nearly inextensional vibrations that can occur in negative -curvature shells 
and the attendant reduced effective stiffness associated with this type of vibration can be 
avoided with proper edge restraints. 
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Percent meridional rise, t 

Figure 7.- Effect of meridional curvature on the minimum frequencies of freely supported 
nearly cylindrical shells for several levels of external hydrostatic pressure. The 
meridional wave number m is 1 unless otherwise noted; the circumferential wave num- 
ber n is given on curves. 


Minimum nondimensional frequency, 


CS 



Nondimensional external pressure, a 

Figure 8 •- Effect of external hydrostatic pressure on the minimum frequencies of several 
freely supported nearly cylindrical shells. The meridional wave number m is 1; 
the circumferential wave number n is given on curves. 


The lowest natural frequency for the class of unstressed freely supported shells 
investigated occurred for t = -3.85 percent with m = 1 and n = 4. (See fig. 2.) Fig- 
ure 9 has a comparison between the minimum frequencies of this shell for each n from 
1 to 15 with freely supported edges and the minimum frequencies of the same shell with 
meridional restraint at the edges (u = 0). From equation (23), inextensional behavior 
occurs for the freely supported shell when m/n = 1/4. The open symbols in figure 9 at 
m = 1, n = 4; m = 2, n = 8; and m = 3, n = 12 locate the inextensional frequencies 
for these particular modes. The inextensional frequencies closely approximate the mini- 
mum for each m -branch of the envelope for the shell with freely supported edges in the 
lower n-range. The in-plane edge restraint had a sizable influence in preventing the 
inextensional vibration and in increasing the frequency for moderate values of n but 
diminished in influence for high and low values of n. The lowest natural frequency 
shifted from m = 1, n = 4 to a more complex mode m = 2, n = 10 and increased by 
an order of magnitude. 

The results for different combinations of edge restraint of the in -plane displace- 
ment u and the slope at the edge w ^ given in table I for r= -3.85 percent show that 
the natural frequencies are highly dependent on the displacement constraint. The slope 
restraint w ^ = 0 has essentially no influence on the natural frequencies. 

Figure 10 is a plot of the minimum frequency envelope as the curvature varies for 
shells with unrestrained in-plane displacements (freely supported) on each edge and with 
restrained in -plane displacements (simply supported) on each edge. The positive - 
curvature shells evidence only a small increase in minimum frequency due to the edge 
restraint, and the effect decreases as curvature increases. The negative -curvature 
shells on the other hand show a large increase in frequency with in-plane edge restraint 
and the dramatic drops in frequency in local regions of negative curvature are no longer 
evident since the in -plane restraint has prevented the development of inextensional 
vibrations. 

For small negative curvature and in-plane restraint there is an initial reduction in 
frequency from that of the corresponding cylinder with a maximum reduction of approxi- 
mately 8 percent occurring at r= -0.7 percent. As the negative curvature increases 
beyond this point, frequency increases and for r = -1.6 percent the shell becomes stiffer 
than the corresponding cylinder. 

Although the analysis and conclusions are restricted to simple shell configurations, 
the behavior of these shells can be used to predict general phenomena in shell structures 
of more complex configurations. For example, it can be expected that local negative - 
curvature areas of more complex shells can develop nearly inextensional modes of vibra- 
tion if the adjoining structure does not provide sufficient in-plane support. Further, 
in-plane compressive meridional stresses can cause severe reductions in the local 
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Figure 10.- Comparison of the effects of meridional curvature on the minimum 
frequencies of freely supported and simply supported nearly cylindrical 
unstressed shells. The number of meridional half waves m is 1 unless 
otherwise noted; the circumferential wave number n is given on curves • 


TABLE I.- EFFECT OF IN- PLANE AND ROTATIONAL EDGE RESTRAINT ON THE 
NATURAL FREQUENCIES OF A NEGATIVE CURVATURE SHELL 
[r = -3.85 percent; s/R = 2; R/h = 50oJ 


n 


ft 



Freely supported 
N^ = v = w = = 0 

(a) 

Clamped with no 
meridional restraint 
N| = v = w= W| = 0 

(a) 

Simply supported 
u = v = w = M£ = 0 

(a) 

Fixed 

u = v = w = w | = 0 
(a) 

i 

0.5278^) 

0.5278( 1) 

0.5709( 1 ) 

0.571l( 1) 

2 

.2429^) 

.2430( 1 ) 

.3663( 1 ) 

.3665( 1 ) 

3 

.0847W 

.0852( 1 ) 

.2532( 1 ) 

.2534( 1 ) 

4 

.00969^) 

,0134( 1 ) 

.1879( 1 ) 

.188l( 1 ) 

5 

.0493^) 

.0502( 1 ) 

.1506( 1 ) 

.150?( 1 ) 

6 

.0782^) 

•0787( 1 ) 

.1307( 1 ) 

.1308( 1 ) 

7 

.0497(2) 

.0532( 2 ) 

.1218( 1) 

.1218( 1 ) 

8 

. 0416( 2 ) 

.0457(2) 

.119o( 1 ) 

.mo* 1 ) 

9 

.0582( 2 > 

.061l(2) 

.1153(2) 

.1154(2) 

10 

.0786( 2 ) 

.0807(2) 

.1118(2) 

.1118(2) 

11 

.085o( 3 ) 

.0893( 3 ) 

.115o(2) 

.115l(2) 

12 

.0949( 3 ) 

.0987(3) 

.1229(2) 

.1229(2) 


a Integer enclosed by parentheses represents meridional wave number. 


stiffness of negative -curvature regions of shells which have only slight in -plane restraint. 
These results suggest that a proper design for a shell structure with negative -curvature 
regions should contain either sufficient in-plane stiffness, such as intermediate rings or 
doubler strips at the boundaries of these regions, or appropriate loadings which provide 
a sufficient positive in-plane stress field to stabilize such regions. 

CONCLUDING REMARKS 

The effect of the interaction between slight meridional curvature and static stress 
due to simple loading conditions on the minimum natural frequencies of nearly cylindrical 
shells is investigated. The effect of various edge restraints is also investigated. 

Axial tension increases the minimum frequencies of negative -curvature shells but 
because of compressive circumferential stresses introduced by this loading in positive- 
curvature shells, the minimum frequencies of these shells are reduced and for sufficiently 
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high tensile loading these shells can become unstable. Axial compression decreases the 
minimum frequencies of negative -curvature shells but increases only slightly the mini- 
mum frequencies of positive -curvature shells. The negative -curvature shells become 
unstable at stress levels well below that of the corresponding cylinder, whereas positive - 
curvature shells offer no increase in stability and become unstable at the same stress 
level as the corresponding cylinder. 

Internal pressure increases the minimum frequencies of both positive- and negative- 
curvature shells and decreases the effect of curvature on the frequency. External pres- 
sure accentuates the destabilizing effect of negative curvature by reducing the minimum 
frequencies, whereas it only slightly reduces the minimum frequencies of positive- 
curvature shells except near instability stress levels. 

In-plane meridional restraint on the boundaries has only a slight effect on the min- 
imum natural frequencies of positive -curvature shells with moderately small positive 
meridional curvature with the effect increasing as the degree of curvature diminishes. 

The in-plane meridional restraint prevents inextensional vibration from developing in 
negative -curvature shells, and as a consequence negative -curvature shells are highly 
sensitive to in-plane restraint, evidencing large increases in frequency above that 
obtained with unrestrained edges. In fact, for a sizable range of curvature the restrained 
negative -curvature shells have higher minimum natural frequencies than the corre- 
sponding cylinder. Rotational restraint has a negligible effect on the minimum frequen- 
cies of both positive- and negative -curvature shells. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 2, 1971. 
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APPENDIX A 

COEFFICIENTS IN EQUATIONS (6), (11), (12), AND (20) 


The coefficients of equation (6) are defined as follows: 
\2 


a 0 = 


- f l( f 7) +f l f 5 f 10 
*4*8 


a 2 


*3(2* 2*7 - *3*5) " 2flf6f7 - (*7) 2 +fl*5*9 + *10^5 +fl*4 - (*2) 2 " 


a 4 = 


H = 



+ *4*10 


*4*8 


-(* 6 ) 2 + f 8^5 +*1*4 - (*2) 2_ 

+ *4*9 


*4*8 


J 


(Al) 


where 


'i - •( i r i£ ) n2 - " 2 


, fl + Ml 
*2 = - n 


f 3 = fr + " 




< 5 =- y » 2 - 2 


f 6 = Xjjt 2 " 


*7 = ~(l + M f-V - A3 


12 
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APPENDIX A - Continued 


f -* 2 

f 8"l2 

f X 2 n 2 N 

f 9 6 " B 


f 10 = 


xV 

12 




N/3 n e 

+ ^n 2 -£1‘ 


The coefficients Uj, Vj, and Wj of equations (11) are as follows: 
For cases 1 and 2, Uj coefficients Uj to U 4 are defined as 


u i “ i 


4 


Re(h]J + Re^hgj 
Re^hj^ — Re^hg^ 


. , at, . B£ 1 


, at . flf 1 
cosh — ^ sin + i 


R* R 


2 2 


I m ( h i) + Im ( h 3 ^ 
S h l) - Ln ( h 3 ) 


. , ai; fit 

sin h cos ~ 

R* R Z 


cosh cos 

R 2 R 2 


U 2-"* 


_ Im ( h l) +Im ( h 3 )] sinh f| sin ^ + 2K hl ) +Re ( h3 )] stoh p cos ^ 

l[ Re ( h l) - Re ( h 3 )] cosh p cos p 


Im ( h l) “ Im ( h 3 ) 


, at . ft 
cosh sin i-tt + 
R 2 R 2 


U 3 = - Re(hg) 


. , (X% . j3f 1 

Slnh p sin ^4 


Im(hi) - Im(hg) 


sinh ^L| cos ^ 
R 2 R 2 


4 


Re/hj'i + Roj'h^ cosh sin ^ + ~ Im^hj) + 


si c„« a 


cosh -5- cos 


R^ 


u 4 = -| 


Im(h]) - Im(h 3 ) 
Im^h ^ + Im^hgj 


sinh sin ^4 + i 


R z R' 


2 2 


. ar| . j3f 1 
cosh — ^ sin 1-2- + 4 

R 2 R 2 2 


Re^hjj - Re^hg^ 
Re^hj) + Re(hg) 


sinh cos ^ 
R 2 R 2 


, af /3£ 
cosh — 71 cos 1-2- 

R 2 R 2 


For case 1 , the coefficients U5 and Ug are 

U 5 = l(h 5 + h 6 )sinh + i(h 5 - h 6 )cosh 


(A 2 a) 


(A 2 b) 


(A 2 c) 


(A 2 d) 


(A 2 e) 
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APPENDIX A — Continued 


U 6 = J( h 5 - he)sinh + I(h 5 + hgjcosh ^ 


(A2f) 


and for case 2, these coefficients are 


U 5 = Re(h 5 )sin ^ + Im(h 5 )cos ^ 

Ue = -Im(h 5 )sin + Re(h 5 )cos 

For both cases 1 and 2, the coefficients U7 and Ug are 
U7 = Re^jsin ^ + 101(117) cos ^ 

Ug = -Im(h 7 )sin ^ + Re(h 7 )cos ^ 


(A2g) 

(A2h) 


(A2i) 

(A2j) 


For the Vj coefficients, replace hj with gj in equations (A2) for Uj. 
For coefficients Wj, the following apply: 

For both cases 1 and 2, 


W 1 = sinh Si sin ^4 
R2 r2 

(A3 a) 

Wo = sinh SA cos — 
R 2 R 2 

(A3b) 

W3 = cosh SA S i n QA 
R 2 R 2 

(A3c) 

W4 = cosh cos ~ 

R 2 R 2 

(A3d) 

Wc = sinh ^4 
b R 2 

(A3e) 

Wg = cosh ^4 
R 2 

(A3f) 
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APPENDIX A - Continued 


For case 2, 


W 5 = S i„^ 

W 6 = cos 2± 
R 2 


For both cases 1 and 2, 


W 7 = sin ~ 

R 2 

W 8 = cos ^ 

Elements of boundary matrices of equation (12) are as follows: 
[n] matrices: 




s 


N 2j = 


nUi + |1 +^)V 4 ’ -^W*’ 


i ’ r l2/ v j " IT J 


?=s 
* 2 


N Si = <j^<2 - |i)V,’ + ^W,'" + - „) - 


N 4j = ^junVj + Wj" - jun 2 W^ 


f=£ 
* 2 


N 5r Uj’ - fxnVj + ^ + |i)w, 


N 6j = 


nU i + l 1+ ll) v i , -4 lw i’ 


t=_s 
* 2 


(A3g) 

(A3h) 

(A3i) 

(A3j) 


(A4a) 

(A4b) 

(A4c) 

(A4d) 

(A4e) 

(A4f) 
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APPENDIX A - Continued 


N 7j = 




. 2 2 
\ n*t* 

“I2" (2 


Nt 

& --r 





(MnVj ♦ W/’ 


- jitn 2 Wj 



where the indices j = 1, 2, . . .,8. 
|V] matrices: 


Y ii = ( D i) . 

Y 2J=( V i) a 

V 'H 

M w l s 

*~2 

Y 4j = ( Wj ’) 

«-! 

y « = ( u i . 

Y 6J - (Vj) 

* "2 

Y 8j = (Wj*) 

£=-S 
5 2 

where the indices j = 1, 2, . . .,8. 


(A4g) 

(A4h) 

(A5a) 

(A5b) 

(A5c) 

(A5d) 

(A5e) 

(A5f) 

(A5g) 

(A5h) 
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APPENDIX A - Concluded 
The coefficients of equation (20) are defined as 


A 0 = b n 


b 22 b 33 " ( b 23) 2 “ b 12( b 12 b 33 " b 13 b 23) +b 13( b 12 b 23 " b 13 b 22) 


2 2 5 

Al = b22t>33 " ( b 23) +b 33 b ll " ( b 13) +b ll b 22 “ ( b 12) 


(A6) 


a 2 = b ll + b 22 + b 33 


where 


b n = fee) + LiJi„2 


b!2' i T ii! f Sn 


b 13 = - 


imrRf R 


\ R 




+ Ml 


b 22 = 




b 23 - - n ( 1 + M £-) - T2 n 
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APPENDIX B 


AXISYMMETRIC VIBRATIONS 


For the particular case of axisymmetric vibrations (n = 0), the circumferential 
equilibrium equation, equation (lb), uncouples from the remaining of equations (1); hence, 
the torsional frequency is independent of u and w. This equation, written in terms of 
the nondimensional meridional length £, becomes 



v ^-f v ,tt = o 


(Bl) 


Since R^ does not occur in this equation, the torsional frequency will be independent of 
the meridional curvature. 

The general solution of equation (Bl) is 


v = 


r 



v. 





j 


(B2) 


so that for the boundary conditions v = 0 the torsional vibration frequencies are given 
by 


Os 


R 





= k7T 


(k = 1, 3, 5, . . .) 


The minimum torsional axisymmetric vibration frequency is thus given by 


7tR 


fyor 




sion 


(B3) 


The two equilibrium equations (la) and (lc) of equations (1) reduce to 


+ u ,tt = o 

+ + ’ 2,1 




^2 ' R^R + j^ W “ B + 13 W ,tt = 0 


(B4) 
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APPENDIX B - Concluded 


These equations are handled in the same manner as were equations (1). The char- 
acteristic roots of equations (B4) are determined from a sixth degree equation formed by 
equating the determinant of the matrix which results after deleting the second row and 
column of the coefficient matrix of equation (12) to zero. 

With these modifications, the terms gj vanish as they should, since there is no 
longer any interdependency between w and v. Since only six roots are present, the 
sum of linear solutions in equations (9) range over six rather than eight terms. With the 
dependence on v removed, the boundary conditions in equations (A4) and (A5) associated 
with v must be deleted, that is, N 2 j, Ngj, Y 2 j, and Ygj, leaving three boundary con- 
ditions on each edge. The solution of the axisymmetric frequencies is determined in the 
same manner as is indicated in the body of the text. Since the shell of interest in this 
report has its torsional frequency as the minimum frequency for n = 0 and is thus inde- 
pendent of r, no calculations for n = 0 have been presented. 
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